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Abstract 



We propose a phenomenological approach to quantum liquids of parti- 
cles obeying generalized statistics of a fermionic type, in the spirit of the 
Landau Fermi liquid theory. The approach is developed for fractional exclu- 
sion statistics. We discuss both equilibrium (specific heat, compressibility, 
and Pauli spin susceptibility) and nonequilibrium (current and thermal con- 
ductivities, thermopower) properties. Low temperature quantities have the 
same temperature dependences as for the Fermi liquid, with the coefficients 
depending on the statistics parameter. The novel quantum liquids provide 
explicit realization of systems with a non-Fermi liquid Lorentz ratio in two 
and more dimensions. Consistency of the theory is verified by deriving the 
compressibility and /-sum rules. 
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I. INTRODUCTION 



There has been increasing interest in generalized statistics for identical particles due 
to possible applications to fractional quantum Hall effect and high-T c superconductivity 
|TJfJ]. Fractional exclusion statistics 0, for which a single-state distribution function can be 
defined in the thermodynamic limit [f|-f| , has appeared in a system of anyons in the lowest 
Landau level J7|, in Heisenberg || and Schrodinger |J quantization of identical particles 
(as fractional statistics in one dimension), in Calogero-Sutherland models |I(J and received 
much attention because of its relevance to the fractional quantum Hall effect, both to the 
bulk [[yj and edge |TJj physics. 

It has been realized that the idea of a generalized exclusion principle || is flexible enough 
to be able to accommodate some other statistical distributions |I3|, for instance those for 



Gentile parafermions [13]. All these distributions belong to a class — we refer to the relevant 
statistics as local — for which the average number of particles in a given quantum state i 
depends only on the Gibbs factor Xj = e^" £i >' kBT for the same quantum state, rij = n(xj). 
The function n(xi) is considered to be quite arbitrary up to natural physical constraints such 
as the existence of the Boltzmann limit for any statistical distribution for an ideal quantum 
gas when the mean occupation numbers are small, — > as Xj — > jjTJJ. If one specifies 
the equation for the distribution function to be 

then the Boltzmann limit condition implies ^{n) — > 1 as n — > 0. 

Low temperature thermodynamic quantities were calculated for ideal gases of particles 
obeying exclusion and other local statistics jlBHTH] • To discuss possibilities to detect gener- 
alized statistics in experiments, one normally has to have a theory of interacting particles. 
In this paper we suggest a phenomenological way to treat interaction in systems of particles 
obeying local statistics, in the spirit of Landau's approach to the Fermi liquid |1| (see also 
p0|). We restrict our attention to statistics of a fermionic type, whose zero temperature 
distribution functions are step functions (with the height of the step generally different from 
unity) . 

In Sect. we formulate the notion of quantum liquids of particles obeying generalized 
statistics and discuss low temperature equilibrium properties for spinless particles. Explicit 
formulas are worked out for exclusion statistics. In Sect. 0I| we study non-equilibrium 
quantities in a simpler case of noninteracting exclusion statistics particles in the relaxation 
time approximation. We show that the Wiedemann-Franz law holds, with a non-Fermi 
Lorentz ratio depending on the statistics parameter, and discuss relation to the calculation 
of the Lorentz ratio for Luttinger liquids |^TJ. In Sect. [TV] we introduce the transport (kinetic) 
equation for the liquid and show that the Lorentz ratio is not modified by the interaction. 
To verify the theory, we check in Sect. [V] that the sum rules for the density response function 
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are respected. In Sect. [VI] we discuss how spin of the particles can be included. We conclude 
in Sect. |V11| with remarks on generic local statistics, non-linear dispersion of quasiparticles 
and possible applications to compressible quantum Hall states. 



II. LIQUID OF SPINLESS PARTICLES 

For an ideal gas of particles obeying fermionic local statistics, by definition, the zero 
temperature distribution function is constant (not necessarily equal to unity) for p < p 
(inside a pseudo-Fermi sphere) and vanishes for p > p . In generalizing Landau's approach 
|H| to such statistics, we assume that when interaction between particles is adiabatically 
switched on, there exists a liquid state which can be described, at low temperatures, in 
terms of elementary excitations, quasiparticles ("dressed" particles), which inherit quantum 
statistical properties of particles of the ideal gas. 

The latter means that the entropy of the liquid, in terms of the quasiparticle distribution 
function n p (we consider spinless particles), is given by the same combinatorial expression 
as for a corresponding gas. For exclusion statistics this implies 



S = J2{[l + (l-9)n p ] Hl + (l-g)n p ] 
v 

-(1 - gn p ) ln(l - gn p ) - n p In n p | . (2) 

The statistics parameter g is assumed to be positive to avoid the Bose limit; g = 1 corre- 
sponds to Fermi statistics. 

The energy of the system is a functional of the distribution of quasiparticles, E = E[n p ]. 
Variation of the energy 

SE = ^2e p 5n p + ^2f pp/ 5n p 5n p/ (3) 
p pp' 

defines the energy of a quasiparticle e p (which is a functional of the distribution function, 
£ P = £p[%>]) as well as the quasiparticle interaction functional / pp / = f pp >[n p ]. 

Varying the entropy at constant total energy and constant particle number, with the use 
of (|3p, yields an equation for the distribution function 

D (p-e p [n p ])/k B T 



e (p-e p \n p \)/k B T _ ^ 



(I - gn p )a[l + (I - g)n P Y-v 
At zero temperature the distribution function is the same as for an ideal gas, 

n p = { l l g \ P \ < Po . (5) 
p I |p| > po v ' 

This can be justified with the use of Eq. (^) and the zero temperature dispersion of quasi- 
particles (H). 
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The number of quasiparticles is equal to the number of particles and hence the relation 
between the pseudo Fermi momentum and the density is the same as for an ideal gas, 

r d d P V 1 

J\ P \< P0 (2ith) d g ' 1 ' 

With d d p = p d ~ l dpdQ p , where dfl p stands for integration over angles of p so that J dfl p = 
Sd-i, and Sd = 2ii d / 2 /T(hd) is the surface of a d-dimensional sphere of unit radius (So = 2), 
Eq. (H) implies a relation between the pseudo Fermi momentum and the particle number 
density p = N/V 

Po \ d pd , , 

\m =9 s d - 1 - (7) 

Note that in microscopic theory, Eq. (^) has be proved, as a generalization of the Luttinger 
theorem [J22|J . 

For T = the quasiparticle energy functional becomes a function since the distribution 
function is fixed. We assume the standard linear dispersion for quasiparticles near the pseudo 
Fermi energy 

s p = ii + v (p-p ) + • • • (8) 

where vo = po/m*, and m* is an effective mass. The quasiparticle interaction functional 
also becomes a function. We will only need to know the quasiparticle interaction function 
on the pseudo Fermi surface, f PoP > - For a spherically symmetric case, it depends only on the 
angle 6 between p and p'\ f PoP > = /(#), with the function f(0) even, due to the symmetry 
fpp 1 = fp'p- 

A relation between the effective mass and the quasiparticle interaction function follows 
from the Galilean invariance. The latter implies J2 P P n p — J2 P m (d£ P [n p }/dp)n p or, upon 



varying w.r.t. n p [[19 



^ = ^-£W^. (9) 

m op dp' 

Multiplying (||) by p, choosing p to belong to the pseudo Fermi surface and using the zero 
temperature distribution function with dn p //dp' = —(p ! /p')^S(p' —po), we obtain 

777* I 

— = i + _„(£„) /(0) cose . (io) 

m g 

The bar denotes the average over angles in spherical coordinates in a d- dimensional p'-space, 
with the polar axis along p and polar angle 9, A(9) = J A(9) (dVL p i / Sd-i)- We have also 
introduced the density of states on the pseudo Fermi surface 
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One can introduce a dimensionless quasiparticle interaction function F{6) = u{eo)f{6), 
and define Landau parameters F and Fi by 

1 



F(9)=F Q , F(6)cos6 = -F 1 . (12) 

In one dimension, where the polar angle takes two values, 9 = 0, 7r, F and F\ are given by 
Fq = |[F(0) +F(tt)] and = |[F(0) — F(tt)]. In two dimensions they are coefficients of an 
expansion F{6) = Fi cos(/6 l ), and in three dimensions they correspond to the standard 
expansion in Legendre polynomials, F{6) = Y^iZo FiPi(cos 9) . 
In the above notation, (|T0|) reads 

- = l + i§. (13) 

For calculating the specific heat one can use the dispersion law (§) for zero temperature. 
Finite temperature corrections to the dispersion (from the second term in the r.h.s. of (El)) 
give terms of higher order in the temperature and can be neglected. The specific heat is 
thus obtained by the change m — > m* in the expression for an ideal gas given in Ref. [17] , 



C v =^v(s )k%T. (14) 

To calculate the compressibility k at zero temperature, defined as 1/k = Np(d/j l /dN)v, 
we note that in varying iV there are two sources for changing the pseudo-Fermi energy, 
Ep[n p \T=o}\p=p , due to the changes of the distribution of quasiparticles and the pseudo- 
Fermi momentum. In this way we obtain 1/k = Np(g + Fq)/v(eq), or in terms of the first 
sound velocity, 

^ 2 = T^_fl + -V (15) 
a mm* \ q I 



Po 



/m\fd in the limit of vanishing interaction). 



III. TRANSPORT PROPERTIES IN THE RELAXATION TIME 

APPROXIMATION 



Four transport coefficients relate electrical and thermal currents to changes in electro- 
chemical potential and the temperature (we consider an isotropic liquid) [|3|] : 

3 = C 11 (E-&)+ £ 12 (-VT) , 
f = C 21 (E - ^) + C 22 (-VT) . (16) 

We start with discussion of a simpler case of the transport equation for noninteracting 
exclusion statistics particles in the relaxation time approximation 
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allowing for dependence of the relaxation time r(e) on the energy. Here n° is an equilibrium 
(uniform) distribution. We do not specify the dispersion law of particles, only assuming that 
it is isotropic, e p = e(p); v(e) = de/dp is the particle velocity. The carrier charge e may be 
different from the charge of an electron and may have either sign. 

Consider the distribution function with slowly varying chemical potential and tempera- 
ture n = n (e^ J '^~ £p ^ kBT ^ r n . Then for the stationary solution to first order in gradients of 
the electrochemical potential and temperature we obtain from ( |17D 



n = n° — t(e 



[v.{eE-V^ + ^{v.VT)Y n ° 



J de 



With this distribution function the transport coefficients in fll6| ) read C 11 = Co, C 12 
-hC 21 = 4=;£i, and C 22 = -hrC 2 , where 



dn° 



rao (in 

-jf de (e - fiy a F (e) — (18) 
and 

a F (e) = e 2 r(e)^25(e - e p )v 2 (e) 
p 

= e 2 T(e)u(e)v 2 (e) (19) 

(v(e) is the density of states) is the energy-dependent conductivity for a Fermi gas with the 
Fermi energy e - The symmetry between C 12 and C 21 manifests Onsager's principle. 

At low temperatures integrals of the form ([Tjf), involving the distribution function for 



r°° 1 rn 

/ den„(e)Q(e) = - deQ(e) 
Jo q Jo 



exclusion statistics, can be evaluated using the expansion |L7 

1 /"A* 

g Jo 

+^Q'^){k B T) 2 + --- (20) 

valid for functions Q{e) which vary on the scale Eq. In this way we obtain the leading terms 
of the transport coefficients 

1 n 2 k 2 T 

\ / e=eo 

The low temperature conductivity a = C 11 = ^<7f{eq) is independent of the temperature. 
The thermal conductivity relates the thermal current to the temperature gradient j = 
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—KVT provided j = 0, thus yielding K = C 22 — C 12 C 21 /C U . From (pIQ we observe that 
the ratio of the thermal conductivity to the electrical conductivity is proportional to the 
temperature thus respecting the Wiedemann-Franz law but with the Lorentz number 

= gL F , (22) 

differed by the factor g from its value Lp for a Fermi liquid. 

This observation tells us that the difference of the Lorentz number from its value for the 
Fermi liquid may signal a kinematic effect (change of statistics). 

The thermopower Q = C 12 / C is 




Q = ~9 



7r 2 k 2 B T (d\na F (e) 



3e \ de 



(23) 



£=£0 



For a free gas, assuming that r(e) oc e x , one obtains ap(e) oc e x+D / 2 and reduces to 
Q = -g(X + d/2)(n 2 k B /3e)(k B T/e ). 

Relation to a Luttinger liquid. — A Luttinger liquid is characterized by the velocities 



of charge and current excitations, and v j, and the sound velocity v s , related by v 2 s = v^vj. 
As two independent parameters, one can choose v s and a parameter /C such that vj = K.v s 
and vn = v s /1C. 

By comparison of the low temperature specific heat for a (spinless) Luttinger liquid 
C/V = 7r 2 kBT/3hv s [^] with that of an ideal gas of exclusion statistics particles, (0) for 



Fq — Fi — 0, we observe that the thermodynamic quantities for the two systems are identical 
if 

g = 1//C (24) 



(see also an argument by bosonization [ 2S| ) . 

Our calculation of the Lorentz ratio show that the the correspondence between a Lut- 
tinger liquid and an ideal gas of exclusion statistics particles can be extended to non- 
equilibrium properties. Namely, by comparison (^) with the result for the Lorentz ratio for 
a Luttinger liquid L = Lp/JC we see that the two expressions coincide under the above 
identification 



IV. TRANSPORT EQUATION FOR THE LIQUID 



Now we turn to the liquid. Consider distributions corresponding to small departures from 
the equilibrium quasiparticle distribution function, n p (r,t) = np(e°) + 5n p (r,t), where the 
superscript labels equilibrium quantities. For short range interactions the corresponding 
change in energy of a quasiparticle can be written in a local form (i.e. depending only on 
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5n p {r, t) at the same point): Se p (r, t) = Y,p> fpp 1 <5%,/(r, t) [pOj] . If n° is expressed in terms 
of the actual energy of quasiparticles, e = e° + 5e, the perturbed distribution function reads 

n p (r, t) = n° p (e p ) + 5n p (r, t) , (25) 

where 5n = 5n — 5e (dn / de) is the quantity (rather than 8n) determining the charge and 
thermal currents, j = —eJ2 p vSn and j Q = —eJ2 P (^ — fi)v8n. The transport equation 

d5n de°d8n dSsdn 
dt dp dr dr dp 

has the same form as for the Fermi liquid, except for the collision integral I(n), where factors 
(1 — n) accounting for the Pauli exclusion principle should be modified. 

The form of the collision integral for ideal gases with arbitrary local statistics for identical 
particles was discussed by the author Hl5| . It was argued in Ref. [fL5|] that in constructing 



collision integrals one should replace factors accounting for the Pauli exclusion principle as 
follows: 

(l-n)->^(n), (27) 

where the function T{n) arises in the equation for the equilibrium distribution function (|]) 
For fractional exclusion statistics the function T{n) can be read off from Eq. (f|): T{n) = 
(1 — gn) 9 [l + (1 — g)n} 1 ~ g . The latter form was also found in discussions of the Boltzmann 



equation for exclusion statistics |26 



Using the change ( p7[) we can write down the collision integral for the case of the residual 
interaction due to an elastic scattering by impurities: 

I{n p ) = -r- Y w pp'$( £ p - s p') {npfinp 1 ) ~ n pl F(n p )} . 
p' 

Here ^■w pp >8(e p — e^) is the transition probability for the scattering of a quasiparticle from 
state p to state p' . We assumed that the impurity atoms are randomly distributed and form 
a dilute gas so that they scatter independently; the detailed balancing principle = w p > p 
was also used. 



Using ( 25|) we linearize the collision integral to obtain 



27T 

I{Sn) = -^-^2w pp/ 5(e p - e p >) [5n p - 5h p/ ] . (28) 
p' 

The transport coefficients can be calculated using the linearized transport equation which 
to first order in gradients of the temperature and chemical potential reads 



eE + ^VTy v ^ = I(5n), (29) 
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where v = de°/dp. 

For an elastic scattering (in particular, for collision integrals of the form fl28|) when w pp > 
depends only on the angle between p and p') , collisions move particles on the constant energy 
surface, and one can seek the solution of the transport equation in the form (cf. relevant 
discussion for a Fermi liquid P7| ) 5n = — (dn° / de)[eE + (p, — e)(VT/T)] • l(p), where l(p) 
satisfies the equation 1(1) = —v. Note that the Landau parameters do not enter explicitly 
(they are absorbed in 5n). As a result, we obtain the same form of the transport coefficients 
C ab as in the previous section for an ideal gas, with the energy dependent conductivity ( |TT?D 
replaced by a pie) = e 2 J2 P $( £ ~ E P ) I 4 v - 



Then calculating the low temperature expansion as in Sect. [[I]] yields the result (pT|). 
with the above crp(e), and leads to the previously found Lorentz ratio (|22|) . 

The fact that the interaction does not change the Lorentz number (for a residual con- 
ductivity due to scattering by impurities) suggests that this number may be a good quantity 
to check in searching for generalized statistics in experiments. 

V. LINEAR RESPONSE AND SUM RULES 

To derive linear response functions from the transport equation, a harmonic perturbation 
of the density by a scalar potential <p(q,u) is applied corresponding to an external force 
F = — iqLp(q, ui)e % ^ qr ~ wt ' ) + c. c. acting on a liquid ||20|| . This modifies the transport equation 
(p6|), written for the Fourier components, to 

(u - qv)5n p (q, u) + qv — ^ f PP i5n p ,(q, u) 

p p 1 

+ qv (p(q, ^)^- = H(8n p (q, lo)) . (30) 

By solving this equation for small (p(q,u), one can find the density response function 
x(q,uj) = (p(q,u))/ip(q,w), where (p(q,u>)) = Y^ P <5n p (q, uS) is the (average) induced density 
fluctuation. 

Consider the collisionless regime where one can neglect the collision integral. In the static 
limit u qv we find from (|30"D 5n p (q, 0) = cp(q, 0)(dn° /de p )/(l + F /g). This leads to the 
response function 

x(«>o) = -irF- (3D 

By comparison with the expression for the first sound velocity (|T5|), we see that Eq. (|3T| ) is 
consistent with the compressibility sum rule liniq^o x(<?) 0) = —N/mv 2 s p0| . 

We note that the transport equation for the liquid is itself valid only for perturbations 
5n varying on a macroscopic scale implying that eq. (BUI) holds only for q small enough. 
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In the quasihomogeneous limit, qvo <C co, in calculating 5n p (q,uj) in (|30|) one has to keep 
terms up to second order in qvo/u, resulting in x(q,c<j) = Nq 2 /mu 2 , which manifests the 
/-sum rule as uj — ► oo J2(J . 

In the hydrodynamic regime, where u is much less than the collision frequency and 
q£ <C 1, with I the mean free path, collisions dominate. In this case, all the deviations 
from the equilibrium distributions are damped out by the collisions, except for those of the 
form (near the pseudo- Fermi surface) 5n^ (constant) and Sn^ cos cp, where cp is the angle 
between k and p. The latter contribute into the particle number and momentum fluxes 
respectively and for them the collision integral vanishes due to the particle number and 
momentum conservation. 

Integrating Eq. (j30|) over momenta, with q as the polar axis, with the weights 1 and p 
respectively, one obtains two coupled equations for Sn^ and Sn^ solution of which results 
in x(Qj u ) = {N/m)/(u 2 /q 2 — v 2 ) leading to the same static and homogeneous limits as 
in the collisionless case. We thus conclude that our phenomeno logical theory respects the 
compressibility and /-sum rules both in the collisionless and in the hydrodynamic regimes. 

VI. LIQUID OF SPINNING PARTICLES. SPIN SUSCEPTIBILITY 

Consider spinning particles, of spin 1/2 for defmiteness, with a = ~, — ~ labeling the spin 
projections (we also use a =|, j in subscripts). 

There is an ambiguity with assigning statistics to spinning particles. In general, one can 
introduce a statistics matrix g aa > (four parameters) determining how particles with different 
spin projections affect each other. There are two special cases, g aa i = gS aa i and g aa i = g 
involving a single statistics parameter. In this paper we discuss the case g aa i = g5 aa > which 
means that particles of different spin projections are not connected statistically (like ordinary 
fermions which are recovered for g = 1). 

Now the energy of a quasiparticle and the Landau interaction function acquire additional 
spin indices. All the formulas of the previous sections are valid upon replacing the quasi- 
particle interaction function by its symmetric combination f{6) — > f s (0) = |[/tt($) + /u(^)] 
(and similarly for the associated Landau parameters) provided, in addition, v(e ) is as twice 
as given by ([□]), due to the two spin polarizations. The latter definition of v{sq) is assumed 
in the rest of this section. 

In the presence of the magnetic field, the variation of the energy of a quasiparticle reads 
(assuming the gyromagnetic ratio two) 

fcpa = ~2fI B aB + fpa,p>a> • ( 32 ) 

p'ot' 

With 5npa = (dn pa /de pa ) 5Epa, the solution of ( |32"D to first order in B at zero tem- 
perature yields modification of the magnetic energy due to the interaction 5e poi = 
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—2\ibolB I '(1 + Ftf/g), where Ffi = ^(F^ >0 — -Fri.o) is an antisymmetric combination of the 
Landau parameters. From the variation of the magnetization 5M = Y^, pa ^^BOtSripa we find 
the magnetic susceptibility 

= A Kfo)_ (33) 

Without including liquid effects (Fq = 0), this expression reduces to the magnetic suscep- 
tibility of an ideal gas of exclusion statistics particles [[18) (which means that the authors 
of Ref. used implicitly the statistics matrix of the form g aa i = g5 aa '). Note that the 



Wilson ratio R w = (iT 2 k] 3 TxV/3{i 2 B C), which is equal to R w = l/(g + F a ) according to (|14D 
and (|33|) , differs from its Fermi gas value (unity) due to both the statistics of particles and 
the interaction. 



VII. CONCLUDING REMARKS 

We have developed a phenomenological approach to a quantum liquid of particles obeying 
fractional exclusion statistics. For generic local statistics of a fermionic type, determined by 
(H), the starting basic equation for the distribution function of quasiparticles (§) should be 
replaced by 

U P = e (v-e P [n p ])/k B T _ ^ 

•F (jip) 

We note that in fact all the formulas for the zero temperature quantities in Sections |T| and 
[VT] as well as the formulas for the linear response (Sect. [V]) are valid in this more general 
case provided g is understood not as the exclusion statistics parameter but just as the height 
of the step determining the distribution function of quasiparticles at zero temperature (|5|). 

The form of the transport equation (|2B[) also holds for generic local statistics. We stress 
that the Lorentz ratio — when it is determined by the scattering by impurities — only contains 
information about the statistics of quasiparticles thus being a good tool for search for new 
statistics. In contrast, the Wilson number encodes information both on interaction and 
statistics and cannot be interpreted in that simple way. 

We have discussed the usual, linear dispersion of quasiparticles near the pseudo-Fermi 
surface. In fact, one can consider more complicated dispersion laws. We comment in this 
respect on a possible application of ideas of generalized quantum liquids to models of two- 
dimensional fermions coupled to a gauge field that were used to describe metallic non-Fermi 
liquid states. 

Specifically, we consider a system of fermions coupled to an abelian Chern-Simons field 
and interacting with a two-body potential v(r) oc Vo/r 11 (or in Fourier components, V(q) = 
Vo/q 2 ~ v ), with l<r]<2(r]=l corresponds to the Coulomb interaction). This model 



11 



is relevant to compressible quantum Hall states with filling factor v = \ |28| as well as to 
high-T c superconductivity fl2~9fl . 

In an attempt to fit the result of the random phase approximation for the above model 
into the quasiparticle picture in the framework of the Landau Fermi liquid theory, a renor- 
malized energy of quasiparticles near the Fermi surface has been introduced, of the form 
28| 



e p - fi oc \p -p F \ 2 sgn(p-p F ). (35) 

2 

The associated specific heat has the form C = ^T^+v . However, the coefficient 7 differs from 
the result obtained in another way, from the free energy of the gauge field (the latter method, 
being gauge invariant, is regarded to be exact) ||30|| . The reason for such inconsistency is 
in that the Landau quasiparticles in this system are ill-defined since the imaginary part of 
the self-energy of quasiparticles is of order of their real part (except for 77 = 1, where the 
damping of quasiparticles is logarithmically small compared to their energy, in which case 
the two ways to calculate the specific heat agree). 

For the quantum liquids discussed in the present paper, if quasiparticles have the disper- 
sion ( [55]) (with pf — > po), the temperature dependence of the specific heat will be the same 
as for a Fermi liquid, with the coefficient 7 depending on the statistics of quasiparticles. 
Based on this, we point out a principal possibility to reconcile the results of calculations 
of the specific heat in the above system in the two ways, by assuming that quasiparticles 
(for i] 7^ 1) form a non-Fermi quantum liquid, corresponding to some generalized statistics. 
In more general terms, our conjecture is that fluctuations of the gauge fields may lead to a 
kinematic effect (changing statistics of excitations from Fermi to some other) for metallic 
states. More detailed discussion of this issue is left for a future publication. 
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